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Abstract — Two-way relaying promises considerable improve- 
ments on spectral efficiency in wireless relay networks. While 
most existing works focus on physical layer approaches to exploit 
its capacity gain, the benefits of two-way relaying on upper 
layers are much less investigated. In this paper, we study the 
cross-layer design and optimization for delay quality-of-service 
(QoS) provisioning in two-way relay systems. Our goal is to 
find the optimal transmission policy to maximize the weighted 
sum throughput of the two users in the physical layer while 
guaranteeing the individual statistical delay-QoS requirement 
for each user in the datalink layer. This statistical delay-QoS 
requirement is characterized by the QoS exponent. By integrating 
the concept of effective capacity, the cross-layer optimization 
problem is equivalent to a weighted sum effective capacity 
maximization problem. We derive the jointly optimal power 
and rate adaptation policies for both three-phase and two- 
phase two-way relay protocols. Numerical results show that the 
proposed adaptive transmission policies can efficiently provide 
QoS guarantees and improve the performance. In addition, the 
throughput gain obtained by the considered three-phase and two- 
phase protocols over direct transmission is significant when the 
delay-QoS requirements are loose, but the gain diminishes at 
tight delay requirements. It is also found that, in the two-phase 
protocol, the relay node should be placed closer to the source 
with more stringent delay requirement. 

Index Terms — Cross-layer optimization, two-way relaying, 
quality-of-service (QoS), delay-bound violation probability, ef- 
fective capacity, resource allocation. 

L Introduction 

The explosive developments of wireless communication 
have brought us into a new era where higher data transmission 
rates and diverse quality-of-service (QoS) provisioning are 
desperately expected. Real-time applications, such as voice 
over IP and video streaming, which are highly delay-sensitive, 
need reliable QoS guarantees. The design merely at the 
physical layer may not ensure the desired QoS requested by 
the service from upper layers. Only through the interaction and 
optimization between different layers can such QoS guarantees 
be fulfilled. This kind of cross-layer approach relaxes the 
layering architecture of the conventional network model and 
brings remarkable performance enhancement, which in turn 
could result in high complexity. Therefore, to develop efficient 
cross-layer methods with small information flows between 
layers is interesting from both theoretical and practical per- 
spectives. 
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Recently, two-way relaying has appeared as an advanced 
relay technique to significantly boost the spectral efficiency 
in wireless networks L2|-[5J. The notion of two-way relaying 
is to apply the principle of network coding at the physical 
layer so that only three or two time slots are needed when a 
pair of nodes exchange information via a relay node, while 
the conventional one-way relaying requires four time slots. 
Most previous efforts on two-way relaying have focused on 
the design and optimization merely in physical layer, such as 
analysis of capacity bounds ||5|, 161, adaptive network-coded 
constellation mapping IT], joint channel coding and network 
coding design ID, ||9l, precoding design with multiple anten- 
nas llTOl . and resource allocation for throughput maximization 
in OFDMA networks |[n]-|[14|. Certainly, it is desirable and 
promising to investigate the cross-layer design and optimiza- 
tion of two-way relay architecture for QoS provisioning. To 
our best knowledge, only two attempts have been made to 
study the cross-layer design for two-way relaying llTSl . llT6l . 
Specifically, IfTSlI and |fT6l characterized the queue stability 
region for infinite backlogs with the XOR and superposition- 
coding based decode-and-forward (DF) protocols, respectively. 
Nevertheless, having stable queue does not provide optimality 
in delay performance. 

Motivated by the needs for bounded delay performance and 
small inter-layer information flows, we adopt the concept of 
QoS exponent and consider the cross-layer optimization of 
two-way relay systems for statistical delay QoS guarantees 
in this work. The QoS exponent is used to characterize the 
statistical delay performance metric, namely, the delay-bound 
violation probability, and is the only requested information 
exchanged between the datalink layer and the physical layer 
l!24l . As a result, our work builds on the information theory 
to capture the performance limits at the physical layer and 
the statistical QoS theory to model the delay performance 
from upper layers. Specifically, we aim to find the opti- 
mal transmission strategies to maximize the weighted sum 
rate while guaranteeing the individual statistical delay QoS 
requirement for each source node. Through integrating the 
theory of effective capacity ll25l . we convert this problem 
into a weighted sum effective capacity maximization problem. 
After that, the optimal power and rate adaptation policies as 
functions of both the network channel state information (CSI) 
and the delay-QoS constraints are developed. 

The main contributions and results of this paper are sum- 
marized as follows: We formulate the cross-layer optimization 
problem for statistical QoS guarantees in two-way relay sys- 
tems as weighted sum effective capacity maximization. This 
problem is shown to be convex. Furthermore, we propose the 
optimal transmission strategies with joint power and rate adap- 
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tation for both three-phase and two-phase two-way relaying 
schemes. Particularly, for the two-phase protocol, the optimal 
channel state partition criterion for successive decoding in the 
multiple-access (MAC) phase is derived. Numerical results 
show that, compared with the conventional two-way direct 
transmission, the considered three-phase and two-phase two- 
way relay protocols can significantly improve the average 
throughput when the statistical delay QoS requirements are 
loose. However, as delay constraints become stringent, the 
performance gain reduces, and eventually all the throughputs 
approach zero. It is also demonstrated that, under the same 
delay-QoS constraint, the three-phase protocol has higher 
weighted sum effective capacity than the two-phase protocol 
in high signal-to-noise ratio (SNR) regime, but has lower 
weighted sum effective capacity than the two-phase protocol 
in low SNR regime. Moreover, we show that the three-phase 
protocol has superiority over the two-phase protocol when the 
relay is extremely adjacent to either of the sources. Meanwhile, 
it is better to place the relay closer to the source with more 
stringent delay requirement for the two-phase protocol. 

The remainder of this paper is organized as follows. Section 
II introduces some preliminaries on statistical QoS guarantees 
and the related work. Section III presents the system model 
and demonstrates our problem formulation. In Section IV and 
Section V, we propose the optimal transmission strategies 
for the three-phase and two-phase transmission, respectively. 
Numerical results are provided in Section VI to verify the 
effectiveness of the proposed policies. Finally, we conclude 
the paper in Section VII. 

II. Background on Statistical QoS Guarantees 
AND Related Work 

A. Preliminaries on Statistical QoS Guarantees 

Due to the time-varying nature of wireless channels, it 
is infeasible to guarantee the hard delay bound for real- 
time traffic. Therefore, statistical QoS metric, in the form of 
the delay-bound violation probability, is commonly used to 
characterize the diverse delay-QoS requirements. 

Based on the large deviation principle, the author in ll24l 
showed that under sufficient conditions, the stationary queue 
length process Q{t) converges in distribution to a random 
variable Q{oo) satisfying that 

lii(Pr{Q(oo) > Qtn}) ^ 



— lim 

Qth->-oo 



Qth 



(1) 



where 6' > is called QoS exponent, denoting the exponential 
decay rate of the distribution, and Qth is the queueing length 
bound. According to the above equation, the probability that 
the steady-state queue length exceeds a certain bound Qth can 
be approximated by 

-SQth 



Pr{Q(oo) > Qth} ~ e- 



(2) 



Similarly, the delay-bound violation probability can be stated 
as 



Pr{D > Dth} « e 



(3) 



the arrival process under a given 6. From (O we can conclude 
that the violation probability for a given delay bound is 
characterized by the QoS exponent 9. Therefore, the dynamics 
of 6 correspond to different delay requirements. Obviously, a 
smaller 6 implies a looser delay QoS constraint, while a larger 
6 means a more strict delay QoS constraint. In particular, when 
6^0, the queueing system can tolerate an arbitrary delay, 
whereas when 6 ^ oo, the system cannot allow any delay. 

Inspired by the theory of effective bandwidth, Wu and Negi 
introduced effective capacity in ll25l . which is defined as the 
maximum constant arrival rate that a given service process 
can support in order to guarantee a statistical QoS requirement 
specified by 9. Analytically, the effective capacity, denoted by 
'^(,(6'), can be given by 

^,(0) = - lim lln(E[e-^^[*l]), (4) 



where S[t\ ~ J2l=i '■^^ time-accumulation of the 

service process, and = 1,2,...} corresponds to the 

discrete-time stationary and ergodic stochastic service process. 
E[-] denotes the expectation. 

Under the assumption that the block fading channel is 
independent and identically distributed (i.i.d) over each time 
frame, the sequence is uncorrected. Then the effective 

capacity can be rewritten as 

1 



<r,(0) = -iln(E[e-''^W]). 



(5) 



where D and Dth denote the queueing delay and delay bound, 
respectively, and (f{9) is known as the effective bandwidth of 



Since the average arrival rate is equal to the average service 
rate when the queue is in steady state, effective capacity can 
also be regarded as the maximum throughput subject to a 
delay-QoS constraint. In particular, as 6' — > 0, the optimal 
effective capacity approaches the ergodic capacity of the 
channel. On the other hand, as 6* — >^ oo, the optimal effective 
capacity is drawing to the zero-outage capacity of the channel. 

B. Related Work on Delay-QoS Provisioning 

Delay-constrained cross-layer optimization has been studied 
extensively in a variety of wireless networks lfT7l - ll23l . lIZTl . 
ll28l . For instance, iflTl focused on the characterization of 
the stability region and throughput optimal control for one- 
way relay systems. The delay minimization problem based on 
power control and relay selection for one-way relay systems 
with multiple antennas was considered in fTSl . Joint power and 
subcarrier allocation for conventional OFDM A networks with 
heterogeneous delay constraints was explored in ||T9l , ||201 . 
Power allocation with statistical delay-QoS provisioning for 
conventional point-to-point, one-way relaying, and multiuser 
systems was studied in ll2Tll . Il22ll and 1231 . respectively. 
Authors in ||27l proposed an optimal scheduling algorithm 
for time division based multiuser systems for statistical delay 
guarantees. Successive decoding order optimization with fixed 
power assignment for MAC channel under statistical delay 
constraints was investigated in |[28|. 

In view of all these existing literature, only the problems 
for unidirectional communication were addressed, while the 
bidirectional nature of the networks has not been fully ex- 
ploited for delay-QoS provisioning. Moreover, the impact of 
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Fig. 1. Cross-layer two-way relay model. 
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(b) Two-phase two-way relaying scheme 

Fig. 2. Two-way relay transmission schemes. 



the cross-layer design and optimization under delay constraints 
in two-way relaying has not been revealed. Therefore, it is of 
great importance and necessity to investigate the two-way relay 
networks for delay-QoS provisioning. 

III. System Overview and Problem Formulation 

A. System Model 

The cross-layer two-way relay system is shown in Fig. [T] 
where two source nodes, A and B, exchange messages via the 
relay node R. Each node operates in a half-duplex manner 
Like in IfTSl , Il22l . we consider that the packets arriving at 
the relay node are forwarded immediately. As illustrated in 
FigH in the datalink layer, two first-in first-out (FIFO) queues 
are implemented at the two sources, which consist of packets 
from upper layer to be transmitted. For QoS provisioning, the 
packets transmitted from one source node to the other are 
subject to the delay constraints, i.e., 6a and 9b- In the physical 
layer, each data packet is divided into frames. Each frame is 
further divided into three or two slots depending on the two- 
way relay protocols to be discussed in Section III-C. 

B. Channel Model 

We consider the scenario in which all nodes have perfect 
channel state information. The channel coefficients of all links 
are assumed to remain unchanged within each time frame but 
vary from one frame to another The instantaneous channel 
coefficient between node i and j is denoted as hij, where 
i,j E {A,B,R} with i ^ j. Here the channel reciprocity 
is assumed, i.e., hij = hj^i, which is valid in time-division 
duplex mode. Without loss of generality, it is assumed that 
the additive noises at all nodes are independent circularly 
symmetric complex Gaussian random variables, each having 
zero mean and unit variance. For notational convenience, 
we further define the instantaneous network channel state 
information as a three-tuple 7 = (71,72,73), where 71 = 
|^A,i?P,72 = |/is,flP,73 = |/iA.sP (as shown in Fig. [T])- 

C. Two-Way Relay Protocols 

Different two-way relay protocols have been studied in 
the literature ||3]-||6l, ||9l. In this paper, we focus on the 
three-phase and two-phase two-way relay protocols with DF 
strategy. Let Pa, Pb, and Pr denote the transmit power of 
nodes A, B, and R, respectively. Let Ra and Rb denote the 
achievable rates from node A to node B and from node B to 
node A, respectively. We further denote P = [Pa, Pb, Pr]^ 
as the transmit power vector, and R — [Ra,Rb]'^ as the 
bidirectional rate pair 



1) Three-Phase Two-Way Relaying: In this protocol, the 
information exchange between A and B is completed in three 
time slots. As shown in Fig. Oa), in the first time slot, source 
node A transmits, while source node B and relay node R 
listen. In the second time slot, node B transmits, while A and 
R listen. In the third time slot, the relay node transmits and 
both A and B listen. Ideally, the time fraction of each slot in 
every transmission frame can be optimized. In this work we 
only consider equal time assignment for simplicity, and so is 
for the two-phase protocol. 

The achievable rate region of the three-phase protocol with 
DF strategy under a given transmit power vector P and 
network CSI 7 is ||5| 

n{P,j) = {{RA,RB)}, (6) 

where 

f \ min {C(7iPa), C(73Pa) + ^(72^^)}, 71 > 73 (7) 
Ra< \ { 

o^hsPA), Otherwise (8) 



f imin{C(72PB),C(73Ps) + C(7iPfl)}, 72 > 73 (9) 
Rb<\ I 

[ -^Cia-iPB), otherwise (10) 

with C(.t) = log2(l + x). It shows that if the channel quality 
of the relay link for one data flow (71 or 72) is worse than 
that of the direct link (73), then direct transmission will be 
triggered for that flow. 

2) Two-Phase Two-Way Relaying: In this scheme, it takes 
two slots to finish one round of information exchange between 
the two source nodes. As shown in Fig. |2b), in the first time 
slot, which is termed as multiple access (MAC) phase, the 
nodes A and B simultaneously transmit signals to the relay 
node R. Due to the half-duplex constraint, there is no direct 
link between nodes A and B. In the second time slot, which 
is known as broadcast (BC) phase, the relay node broadcasts 
signals to both A and B. 

The achievable rate region of the two-phase two-way relay- 
ing with DF strategy under a given transmit power vector P 
and network CSI 7 is ||1-||3 

Tl{P,^)=CMAc{PA,PB,l)f^CBc{PR,l), (H) 
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where Cmac and Cbc are the achievable rate regions of the 
MAC and BC phases, respectively, and can be given by 

CuAciPA.PBn) - {{Ra.Rb) ■■ Ra < ^CijiPA), 

Rb < Ic{^2Pb), Ra + Rb< ^C(7iPa + 12Pb)}, (12) 



CBciPRn) - {iRA,RB) ■■ Ra < ^C{-/2Pr), 

Rb<\c{jiPi,)}. (13) 

Note that both Cmac and Cbc are convex sets, so is their 
intersection. 



D. Problem Formulation 

In this paper, our objective is to find the optimal trans- 
mission policies to maximize the weighted sum rate of the 
two-way relay system while satisfying the individual delay 
requirement at each node. As stated before, the effective 
capacity can be viewed as the maximum throughput under 
the constraint of QoS exponent in steady state. Hence, we can 
formulate an equivalent problem, which is to maximize the 
weighted sum effective capacity for given delay constraints of 
node A and node B, i.e., 6a and 9b- Our resource allocation 
policies are based on cross-layer parameters, specifically, the 
instantaneous network CSI 7 = (71,72,73) and the delay- 
QoS requirements = {9 a, Ob)- Correspondingly, we define 
e = (7,^). Therefore, the problem can be formulated as 
follows. 



PI 



max 

P(£),R(£) 



s.t. 



9a 




9b 


(14) 


E^[Pa]<PI 


(15) 


E^[Pb]<P^ 


(16) 


^■y[PR]<P^ 


(17) 


i^(e)G7^(P(e),7) 


(18) 


Pie) h 0, 


(19) 



where uja,^-^b are the weights assigned to the two users, 
satisfying ua + i^b = 1, and Pa, Pb, Pr are the long-term 
power constraints of node A, B and R, respectively. E^[-] 
emphasizes that the expectation is with regard to 7. P(e) 
and -R(e) denote the power and rate adaptation policies to 
be optimized, which are functions of e. The instantaneous 
rate region TZ{P{e),j) is defined in ^ for the three-phase 
protocol, or in (flU for the two-phase protocol. Note that 
TZ{P{e), 7) is a convex space spanned by the power sets P(e). 

It is proved in 1231 that the weighted sum effective capacity 
is a concave function of the powers in multiuser systems 
with direct transmission. Using the similar method, we can 
prove the weighted sum effective capacity in the two-way 
relay system as given in (fT4l i is also concave of P{e). The 
main reason is that the achievable rate pair Ra and Rb are 
both concave with respect to the power vector P. In addition, 
the power constraints (fT5Tl-(fT7]i and JT^t are affine. Thus, the 




(b) An example of two-phase 
region for 7^'' 



(a) Three-phase region 
Fig. 3. Cross-layer two-way relaying channel state region. 

problem PI is a convex optimization problem, and there exists 
a unique and optimal solution. In the next two sections, we 
shall develop the optimal power and rate adaptation policies of 
the given problem for the three-phase and two-phase protocols, 
respectively. 

IV. Optimal Policy for Three-Phase Two-way 
Relaying 

In this section, we first derive the optimal transmission 
policy for the three-phase two-way relaying subject to general 
delay QoS requirements 9 ~ {9a, 9b), for which TZ{P{e), 7) 
is given in (|6]l. Then we consider the limiting case when both 
6a and 9b approach zero, i.e., the ergodic capacity problem. 

A. Optimal Policy 

We define the Lagrangian of problem PI as 



C{P{e),R{e),X 

UJA 



^ In (E-,[e-'^^^^^^]) - ^ 111 (E^[e-''-«-(^)]) 
6a 9b 

+ XA(PX~^-y[PA{e)]) + Xb{P^ -K-y[PB{e)]) 
+\n{l^-E^[Pn{e)]), (20) 

where A = [Xa, Xb, Xu]'^ are the Lagrange multipliers related 
to the power constraints (fT5Tl-(fT7]i. Then the dual problem of 
PI can be stated as 

P2 : min max C(P(e),R(e),X). 

x>o P(e)yo, ^ ^ ' \ ' / 

R{e)eTZ 

Note that the subgradient method can be used to update A 
toward the optimal A* as follows 



A(') _s('' (P-E^[P(e)] 



(21) 



where the subscript i denotes the iteration index, and is 
the vector of step size designed properly. 

According to the achievable rate region defined in (|6), we 
divide the channel states 7 into four regions as shown in 
Fig. |3ja). In the following, we provide the optimal trans- 
mission policy for each channel region, and the detailed 
derivations are given in Appendix lAl 
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1) Region Ri(7i < 73,72 < 73)-' In this case, the 
achievable rate pair Ra and Rb satisfy (8) and (10), respec- 
tively, which means that the bidirectional links follow direct 
transmission. It is obvious that the optimal rates are exactly 
the capacity bound, given by 



RA{e) = \c{^^PA{e)), 
RB{e)^\c{^^PB{e)). 



(22) 
(23) 



Substituting the above into the Lagrangian function ( l20l i to 
eUminate the rate variables, we can then obtain the closed- 
form expressions of the optimal power allocation as 



PA{e) - 



73 



Pa 



73-^ 



"FST3 



73 



PBie) = 

where a = 3\n2, I3i ^ e {A, B}, and 



73 



(24) 



(25) 
(26) 



Note that 4>i, (j32 are expectations over all the regions, and 
should be updated with A in each iteration of the dual problem. 

2) Region M2(7i > 73,72 < 73)-' In this region, node A 
transmits signals with the help of node R while node B adopts 
direct transmission. According to (7) and (10), the optimal rate 
allocation is given by 

RA{e) = imin{C(7iPA(e)),C(73/'A(e))+C(72Pi?(e))}, 

(27) 



i?B(e) = ic(73PB(e)). 
Then the optimal power allocation can be obtained as 

■pA{e)=[PA{e)]-' 
PB{e) - 



(28) 



^) 73 

_ hl-73)PA{e) 



73 



(29) 



where Pa{^) is the solution of the following equation 



a4>i 



72[l+73^'_A(e)] 

(30) 

We use the simple bisection method to obtain Pa{^) since 
( [3OI 1 is a monotonically decreasing function of Pa{^)- 

3) Region K3(7i < 73, 72 > 73)." Node B transmits signals 
via the assistance of the relay while node A adopts direct 
transmission. This case is similar to that of R2 and we omit 
the results here. 



4) Region 1^4(71 > 73,72 > 73)-' In this region, both 
source nodes need the relay node's help. According to the 
achievable rate pair (7) and (9), the optimal rate should follow 

Ra{€) = \mm{C{jiPA{e)),C{j3PA{e))+C{^2PR{e))}, 

(31) 

Reie) = imin{C(72Pi3(e)),C(73PB(e))+C(7ii^i?(e))}. 

(32) 

The associated optimal power allocation is as follows. 
Define 

A 71(71 - 73) 



72(72 - 73)' 

If r < 1, i.e., 73 < 71 < 72, we have 

■ PAie) ^ [PAie)] + 
Psie) = j^^^ 

Pr{^) 



(33) 



+ (l-r)73PA(e) 

(71 — ri)PA{<^) 

72[l+73fA(e)] ' 



(34) 



where Pa{^) is the solution of 

■[l+7iP^(e)] 

T72^^A(e) 



iOA^ii - ^ - _EA+jl 



TUJBI2 



l + (l-r)73PA(e) 

t\b 

[l + (l-r)73p^(e)] 



cj4>2 [l + (l-T)73P^(e) 

' -^^(71-73) 



72[l+73^'A(e)] 
2 - Aa = 0. (35) 



Note that the bisection method can be used to obtain Pa{^) 
since ( [35] l is a monotonically decreasing function of Pa{^)- 
If r > 1, i.e., 73 < 72 < 7i, we have 



Pb(6)=[Ps(6)] + 

p /,\ _ (72-73)-PB(e) 
^^R{e.) - 7i[l+^3PB(e)] ' 



(36) 



where Psif) is the solution of 



^B72 



l + 72Ps(e)]" 



7i^B(e) 



TUJAll 



r+(T-l)73PB(e) 

tXa 

[r + (T-l)73PB(e)]' 



(701 [r + (T- 1)73^5(6)]' 

,a_4+3 

' Ai?(72 - 73) 



7i[l + 73PB(e)] 
As = 0, (37) 



which can also be obtained by the bisection method as ( [35] ). 

In summary, the optimal transmission policy {P(e), -R(e)} 
for the three-phase two-way relaying is given by (l24t . (l22t . 
(|23] | when 71 < 73,72 < 73, by (|29l), (|27li, (|28l) when 71 > 
73,72 < 73, by (I34I 1. (I3TI 1. ( |32] | when 73 < 71 < 72, and by 
(|36] |. (I3TI 1. ( I32I ) when 73 < 72 < 71. The detailed derivations 
are given in Appendix |A] 
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B. A Special Case 

As reviewed in Section II-A, the dynamics of 9 correspond 
to diverse delay-QoS requirements. In particular, in our two- 
way relay system model, if 9a ^ Ob ^ 0, meaning the 
services of the two nodes are non-real-time, then the weighted 
sum of effective capacity yields the weighted ergodic capacity. 

Letting 9a ^ 9b ^ in (EHi, and we can 

obtain the optimal transmission policy for weighted ergodic 
capacity maximization. For example, in region Mi, the optimal 
powers are given by 



P* (e) = 0, 



to A 



73 



7,3 



(38) 



which have the standard form of water-filling. For the rest 
three regions, similar results can be obtained. 

V. Optimal Policy for Two-Phase Two-way 
Relaying 

In this section, in accord with Tl{P{e),-f) given in ( fTTT ). 
we present the optimal transmission policy for the two-phase 
two-way relaying, as well as the optimal partition criterion 
in the MAC phase. Meanwhile, the problem in limiting case 
when 9a = 9b = i& also analyzed. 

A. Optimal Policy 

The rate constraints in the BC phase of this protocol can be 
equally rewritten as 



1 



1 



> — e 

- 9a 

^B ~ Ob 



-^C(72Pfl(e)) 
-^C(7lPR(e)) 



(39) 
(40) 



By using the Lagrange dual method 11291 , we can involve the 
two rate constraints into the objective function of PI. Then, 
the resulting Lagrangian can be expressed as 



i^(P(e),ii(e),A,/x) 



LOB 



M / „-eARA{e) _ „-^C{r2PM{e)) 

Oa V 



+XB.(P^,-E^[PRie)]), 



[PB{e)] 



(41) 



where /j, ~ [/i^j/is]-^ are the Lagrange multipliers associated 
with the rate constraints on Ra{£) and i?_B(e) in i39[ . ( |40| |. 
A = [Xa , Xb , X]i] ^ are the Lagrange multipliers related to the 
power constraints. As a result, the dual problem can be stated 
as 

P3 : min max ^(P(e),R(e),X,a). 

R{e)eCMAC 



Taking a close look at the above dual function, we find 
that the optimization of the relay power policy PrXc) can be 
decoupled from others. Therefore, the dual function can be 
computed by solving the two subproblems as follows. 



subproblem 1 : 



max 

PA(e)>0, 
PB(e)>0, 
RA{e),RB{e) 



^In (E.,[e- 



iRA{e)] 



iRAie) 



fJ-B 



+XA{Pl-Ey[PAie)]) 

+AB(P^-E^[PB(e)]) (42) 
s.t. (i?A(e),i?B(e)) eCMAc(P4(e),Ps(e),7),(43) 



and 



subproblem 2 : 



max 

PR{e)>0 



E. 



MA 



^(l + 7iP,^(e))- 



l+l2PR{e)y 



XR{PR-E^[PR{e)]), 



(44) 



where = j^,* G {A,B}, is the same as the definition in 
Section IV-A. Then the dual problem P3 can be solved through 
two nested dual searching loops. The inner loop searches 
fi for given A and the outer loop searches A for given fi, 
where the values fi and A can be updated iteratively using the 
subgradient method with guaranteed convergence as 



Ma 



Ma 



-eARAi^) 



.-YC(72PH(e)) 



^1;+') = - 4^ ie-'^M^^ - 6-^^7(7, P„(.)) 



^ a(') - s(') (P - E.^[P(e)] 



, (45) 
, (46) 



(47) 

where the subscript i denotes the iteration index, and zy, 
and s*^*^ are the step sizes designed properly. The overall 
algorithm is specified later In the following, we solve the two 
subproblems respectively. 

1) Solution of Subproblem 1: The first subproblem is 
relevant to Pa{£) and PB(e) in the MAC phase but not PR{e) 
in the BC phase. It is essentially a classical resource allocation 
problem in MAC channels, though the objective function is 
slightly different from the Gaussian MAC ll26l . As shown in 
ll26l . successive decoding is the optimal strategy for resource 
allocation in MAC channels. Motivated by this result, we 
partition the channel states into two regions, R'^ and Rj, for 
which an example is shown in Fig. Ob). In region M.[, the 
relay first decodes the signal from node A, then subtracts this 
decoded signal from the received signal, and then decodes 
the signal from node B. Inversely, in region the relay 
decodes the signal from node B first and then the signal from 
node A. Similar partition method is used in the literature (e.g., 
1271 . ||28| ). In the following, we propose the optimal power 
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and rate adaptation policy for a given channel partition. The 
optimal channel partition method will be derived in the next 
subsection. 

In region M.[, the maximum rates are achieved at 



_ I f 71^^ (e) 
RB{e) = ic(72Ps(e)), 



while in M.^, the maximum rates are achieved at 

RA{e) = ic(7iPA(e)), 



-C 



12Pb{^) 
l + 7iPA(e) 



(48) 
(49) 

(50) 
(51) 



By applying the Karush-Kuhn-Tucker (KKT) conditions 
1, the optimal power adaptation policy for Pa (e) and Pb (e) 
in M']^ must satisfy the following conditions (the derivation is 
provided in Appendix IbTi: 



PA{e) = 
PB{e) = 



WPPl j l+72-PB(e) 
71 



Pa 



l+72PB(e) 
71 



T + 



(52) 



PB{e)r, 



where Pb{^) can be obtained using a numerical searclQ 
through the following equation 



— [l+l2PB{e) 
a2 



Aa72 
Ab7i 



Aa72 / "1 \ '^-**^^ 



Ab7i \7i 



[1 + 72^5(6)] ^^+^-1 = 0,(53) 



with 



Oil 



SXa 



XI -1 ' 



0L2 



SXe 



XI -1 ■ 



l+l2PB{e) 



J A 



+ 



.Ea 



l + 7iPA(e)] ^ p-y (71, 72)^71 ^^72, 



r 1 — 

[l + 72-PB(e)J p-y (71, 72)^71 c^72 



■yel 



72^*5 (e) 



l + 7iPA(e) 



(71, 72)^71^72, 



where S — 2 In 2. Like the three-phase two-way relaying, (p[ 
and 02 are also updated with A in each iteration. 

Using the similar method, the optimal power allocation in 
region can be obtained and the details are omitted. 



The outline of this numerical search is: First, find out the stationary point 
based on the derivative of the function. Then, determine the interval where 
the zero point exists, and search for this point using the bisection method. 



2) Solution of Subproblem 2: The second subproblem is 
only relevant to PR,{e) in the BC phase, which is also a convex 
problem. By applying the KKT conditions, we can obtain the 
optimal power allocation PR.{e) = \_PR{e)\^ , where Puie) 
must satisfy the following equation 

/^A72[l + 72^'/?(e)] 

+^lBll[l + llPR{e)^^ -SXr^O. (54) 

We apply the bisection method to obtain PR{e) since ( l54l i is 
a monotonically decreasing function of PR{e). 

B. Optimal Partition Criterion 

In the above subsection, we have demonstrated the optimal 
transmission policy for given decoding region R' = (R']^,M2) 
in the MAC phase. Here we present the optimal partition 
criterion to determine the decoding order in the MAC phase 
based on the obtained transmission policy. 

To maximize the weighted sum effective capacity in ( l42b . 
the optimal decoding order in the MAC phase should be 
dynamic with respect to different channel state information. 
Similar to ||27l , ||28l . finding such optimal partition can be 
viewed as finding an optimal threshold, 7^'' for 71 or 72'' for 
72. Here 72'' (or 7*'') is a function of both 71 (or 72) and 
the power allocation policy P{e). The following proposition 
is established to find the optimal threshold. 

Proposition 1 ( Optimal Channel Partition Criterion ): 
When 72'' is used to partition M', 7 falls into region Rj if 
72 < tI'', otherwise 7 falls into region Rj, where 72'' must 
satisf\^ 



.f^A 



1+7^"Pb(£)J 



[l + 7iPA(e) 



/3a 



[1- 



7^''-PB(e) 
1+71 Pa (e) 



K. (55) 



l+7fPs(e) 



When 7* is used to partition M', 7 falls into region R'j^ if 
71 > li^, otherwise 7 falls into region Rj, where 7*'' must 
satisfy 



1 



llPA(e) 1 
l+72Pf3(e)J 



_Pa 



1 



lY'PAie)]' 



.fA 



[1 + 



_72Pb(£) 1 ^~ 



[l + l2PB{e)\ 



K. (56) 



l+-tY'PA(e)i 

In both ( |55] | and ( |56] ), K is defined as 



K 



dB{uJA<i>'i -ma) 
Proof: Please see Appendix ICl ■ 
Particularly, 7*'' and 72'' should be well-defined in the parti- 
tion criterion, meaning that they should be positive. However, 
as we know from ( fSSl l and (|56] |, they cannot be both positive 
in some condition. In this case, the obtained positive one is 
chosen for the partition. 

Finally, we describe the overall algorithm in Algorithm 1 
to find the optimal power and rate adaptation policy for the 
two-phase two-way relay protocol. Note that in Algorithm 1, 
for a given channel partition we can obtain the optimal power 

-We use the similar numeiical method as described before. 
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allocation, then the optimal power allocation in turn leads to an 
optimal channel partition. Due to the convexity of the problem, 
the global convergence and optimahty can be guaranteed. 



C. A Special Case 

Similar to the three-phase scheme, when 6a—0b^^,^^ 
can obtain the optimal transmission policy for the two-phase 
two-way relaying without delay requirements (i.e., the ergodic 
capacity). 

Proposition 2: The optimal power allocation policy for the 
two-phase two-way DF relaying for weighted ergodic capacity 
maximization when S^a < is determined by 



£,A £,B - 

SXa ^- '^"^ 



Cb-0 



1 

72 



[0, 

I 2ci ■ 



^ 3 

otherwise 



(57) 



(58) 



(59) 



where ^a = ^A~f^A, ^b = i^B-f^B, and ci = XrJiJ2, C2 = 
A/?(7i+72)-7i72(/^a + M-b)/5, C3 = A_r-(/^a7i+Mb72)/'5. 
The results for £^a > Cb can be easily obtained using the same 
methods. 

Proof: Letting 6a = Ob ^ in (l52]i and ^5^, we have 



71 



L(5(Ab7i-Aa72) 

^.4 72 I M-B71 

l+72PH(e) ^ l+7l-Pj?(e) 



72 J 

5Xb^ 0. 



Through simple calculation, we can get the desired results 
(|57]i, (|5l and ^ in R^. 

For the segmentation of (M']^, Mj) the case of 6'^ = 6'b ^ 
0, by applying the proposed optimal partition criterion, we 
have 

6A{t^B 



K 



Under this condition, dSST l becomes 



l+7iFA(e) 
l + -ifPB{e) 



= K. 



When > S^A, meaning that K > 1, > is well-defined 
for partition. As 9a 0, we can find 



which implies that 7 always falls into Mj. Hence, the proof 
completes. Similar analysis can also be done when < -FT < 1 
by using 7^''. ■ 



Algorithm 1 Finding the optimal transmission policy for two- 
phase two-way relaying 
1: Given lua, i^b- 

2: Initialize PA{e) = PI, PB{e) = PE- 
3: Initialize A. 
4: repeat 

5: Initialize fi for each 7. 
6: repeat 

7: Determine the decoding order in the MAC phase 
according to the optimal partition criterion given in 
Proposition 1. 

8: if the signal from node A is decoded first, i.e., 7 e M.[ 
then 

9: Obtain the two source power allocations Pa (e) and 

Peie) by using 
10: Obtain the relay power allocation PR{e) by using 

dSli. 

11: Obtain the rate adaptation R{e) as (|48]l, ( |49l ). 

12: Update /J, using the subgradient method in ( |45] l, 

13: else 

14: 7 G adopt similar operations as above except 

that the rate adaptation is given in dSOl l. ( BTl i. 

15: end if 

16: untU /J, converges. 

17: Update A using the subgradient method in ( |47] |. 

18: until A converges. 



VI. Numerical Results 

In this section, extensive numerical results are provided 
to illustrate the performance of our proposed cross-layer 
transmission strategies for the two-way relay systems. 

As a benchmark, the conventional two-way channel with 
direct transmission is considered, for which the optimal trans- 
mission policy is obtained from ET\ . Besides, to show the 
advantage of the optimal channel partition in the two-phase 
two-way relaying, the transmission policy using the static 
weight-based channel partition, as introduced in ||261 . is stud- 
ied. Specifically, in the MAC phase of this scheme, the relay 
always first decodes the signals from the source with smaller 
weight, regardless of the CSI variation. The weight-based and 
the proposed CSI-based schemes with fixed power assignment 
are also studied to illustrate the significance of channel-aware 
power adaptation. In these fixed power assignment schemes, 
the instantaneous power of all transmitting nodes is set to be a 
constant, while the transmission rates are adaptive with respect 
to the channel fading. 

In our numerical evaluation, the relay is located in a line 
between the two users. We set the distance between node A 
and B as 2. The A-R distance and the B-R distance are 
denoted as d and 2 — d, respectively, where < d < 2. 
The log-distance path loss model with small-scale Rayleigh 
fading is assumed. Hence, the network channel information 
7i : 72 7 73 follow independent exponential distribution with 
parameter Ai ~ d^ , A2 = {2 — d)'\ and A3 = 2'^, respectively, 
where v denotes the path loss exponent. A typical value 
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User A's effective capacity (bits/s/Hz) 



Fig. 4. Optimal effect ive c ap acity regions for d ifferent scliemes under diverse 
delay constraints with P4 = Pg = 9dB, Pn = 6dB. 



3.5 




,1 , , , , , 1 1_ 

2 5 8 11 14 17 20 23 

Users' long-term power constraints (dB) 



Fig. 5. Weighted sum effective capacity versus long-term power constraints 
under Oa = 9b = 1- 

of ly lies in the range of (2,5), and it is set as 4 in our 
examples. The long-term power constraints of the nodes satisfy 
Pa = Pb — Pr + 3dB. Throughout this section, the weights 
are given by uja = 0.6 and ujb = 0.4 for illustration purpose. 

A. Performance of Symmetric Relay for Two Sources 

In this subsection, we consider the case when relay is in the 
middle of the two source nodes, namely, d — 1. Hence, the 
channels between the sources and the relay are symmetric. 

Firstly, in Fig. H] we plot the optimal effective capacity 
regions under different delay constraints for the three-phase 
and two-phase two-way relaying, as well as the two-way direct 
transmission. We set the long-term power constraints of the 
source nodes as 9dB. It is observed from Fig. |4] that, with 
the help of two-way relay, the effective capacity region is sig- 
nificantly expanded compared with the conventional two-way 




lO"'^ 10"' 10° io' lo'^ io' 

Delay-QoS constraints 6^=0^ (1/bits) 



Fig. 6. Weighted sum effective capacity versus different delay constraints 
with long-term power constraints Pa = Pb = 9dB, Pr = 6dB. 



240 




10"'^ 10"' 10° 10' 10^ 10^ 

Delay-QoS constraints 6^=6^ (1/bits) 



Fig. 7. Effective capacity gain of optimal poUcies over direct transmission 
with long-temi power constraints Pa = Pb = 9dB, Pr = 6dB. 

direct transmission. We can also find that, if one user's delay 
constraint becomes stringent, its effective capacity becomes 
small, so is the effective capacity region. This suggests that 
there is in general a fundamental throughput-delay tradeoff 
associated with the optimal resource allocation. In addition, 
under these given power constraints, the effective capacity 
region of the two-phase protocol is larger than that of the 
three-phase protocol. 

In Fig. |5] we compare the weighted sum effective capacity 
of different schemes under different long-term power con- 
straints, where the delay constraints are set as 9a = Ob = ^■ 
From this figure, we can observe that the two-way relaying 
brings tremendous improvements on the weighted sum effec- 
tive capacity for information exchange between sources. By 
taking a closer look at Fig. |5] it can be found that the power 
adaptation can bring about 10% and 7% improvements on the 
effective capacity for the two-phase and three-phase protocols. 



10 



X 

1 0.8 









2 phass optimal 




2-phase fixed power 


' /I 


3-phase optimal 


3-phase fixed power 


/ / ^ \ 


\ \ 




\ \ 


1 ' / AV 


. \ \ \ 
\ \ \ 


/ // / \ 


\ ^ \ 


/ // / 3— phasG 


\ \ \ 


1 1/ f 




Id/ 






\ \ \ 


A' 


\V,\ 


/// 


\\\ 


/' 


v\ \ 












^ ■ 






/ / 


\ 







0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



A-R distance (d) 

Fig. 8. Impact of relay location on weighted sum effective capacity with 
same delay constraint Q js^ = = \. 



respectively, at all the considered power constraints. For the 
two-phase protocol, about 5% performance improvement can 
be achieved by the proposed CSI-based scheme over the 
weight-based scheme. Moreover, we can also see that the 
weighted sum effective capacity of the two-phase protocol is 
superior to that of the three-phase protocol when the source 
power is below about 18dB, while it is inferior to the three- 
phase protocol when the source power is higher Therefore, the 
three-phase scheme with power adaptation is more appropriate 
for cross-layer two-way relaying under high SNR conditions. 

Fig. |6] shows the weighted sum effective capacity versus 
different delay-QoS constraints d = 6a = Ob, and the 
effective capacity gain of the optimal policies over direct 
transmission are further plotted in Fig. |7] Here, the power 
constraints for the source nodes are fixed as 9dB. As presented 
in Fig. |6] the weighted sum effective capacity generally 
decreases with the increasing 6. We can observe from both 
Fig. |6] and Fig. Q that when the delay constraints are loose, 
the optimal policies for the two-phase and three-phase two- 
way relaying can achieve substantial effective capacity gains 
over the direct transmission. However, the advantages become 
small as the delay constraints go stringent. Particularly, the 
weighted sum effective capacity of all schemes approach to 
zero if 9 is large enough. This is expected as, when the 
delay constraints are very stringent, the system can no longer 
support the transmission due to fading effect of the channel. 
This conclusion is consistent with the theory of delay-limited 
capacity in information theory. Moreover, it is obvious that the 
proposed strategies can efficiently provide the best weighted 
sum effective capacity in two different protocols. Specifically, 
the margins between the optimal policies and the fixed power 
policies for both protocols go larger along with the decrease 
of 6. Besides, it has demonstrated that, under such condition, 
the two-phase protocol is superior to the three-phase protocol, 
though the superiority becomes small for stringent delay 
constraints. 



0.5 




0.1 ' • • • • ' ■ ■ ■ • ' 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
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Fig. 9. Impact of relay location on weighted sum effective capacity with 
delay constraints 8a = 100,6b = 1. 

B. Impact of Relay Location on System Performance 

In this subsection, we will consider the impact of the relay 
location in the two-way relay system. Here, we assume the 
long-term power constraints of the source nodes are 9dB. 

First, we consider the case when source nodes have the 
same delay-QoS constraint that 6*^ = 6'b = 1. As illustrated 
in Fig. [8] the maximal weighted sum effective capacity is 
obtained when the relay is in the middle of the two source 
nodes no matter what transmission strategy is adopted. Mean- 
while, we can observe that our proposed resource allocation 
policies can obviously achieve effective capacity gains in both 
transmission schemes. However, the benefits decrease when 
the relay is close to either of the source nodes. The reason is 
that, as a result of the severe path loss, the channel between 
the relay and the distant source becomes the major limit of 
the transmission in this case. Moreover, we can find that, on 
this condition, when the distance from the relay to the source 
node is less than 0.3, the three-phase protocol outperforms 
the two-phase protocol. Otherwise, the two-phase protocol has 
advantage over the three-phase protocol. 

Next, we study the situation where the two source nodes 
have different delay requirements. Here, we set that 6 a = 100 
and 9b = 1- It is interesting to find from Fig. |9] that, the 
maximal weighted sum effective capacity for the three-phase 
protocol is gained when the relay is in the middle of the two 
sources, while the relay should be closer to the source with 
more stringent delay requirement if the two-phase protocol is 
used. Different from Fig. |8] when the distance between the 
relay and the node with greater 9 is larger than about 1.3, the 
three-phase protocol can get better performance. 

VII. Conclusion 

This paper has studied the cross-layer optimization of two- 
way relaying under statistical delay-QoS constraints. We have 
focused on two main transmission protocols: three-phase trans- 
mission and two-phase transmission. By integrating the theory 
of effective capacity, the optimization problem for weighted 
sum throughput maximization in physical layer and delay 
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provisioning in datalink layer was modeled as a long-term 
weighted sum effective capacity maximization problem. Then, 
optimal transmission policy was proposed for each protocol. 

A few important conclusions have been made through 
extensive numerical results. Firstly, our proposed two-way 
relaying policies can efficiently provide delay-QoS guarantees, 
while there exists a tradeoff between the throughput gain and 
the delay-QoS provisioning. Secondly, the proposed policies 
significantly improve the system performance compared with 
the schemes without power adaptation. Especially, for the two- 
phase protocol, the proposed CSI-based method for succes- 
sive decoding in the MAC phase has 5-10% performance 
improvements compared with the weight-based successive 
decoding. Thirdly, when the relay is located in the middle of 
the transmission, in terms of weighted sum effective capacity, 
the three-phase procotol outperforms the two-phase protocol 
in high SNR regime and is inferior to the two-phase protocol 
in low SNR regime. Last but not least, it is better to place the 
relay closer to the source with more stringent delay constraint 
for the two-phase protocol. 

This work has concentrated on DF relaying with full 
channel information in two-way relay systems. It can be 
further extended to the heterogeneous networks consisting of 
delay-constrained and non-delay-constrained traffics. Alterna- 
tive relaying protocols and transmission strategies with time 
adaptation or partial channel state information are also possible 
for future research. 

Appendix A 

Derivation of Optimal Power Adaptation for 
Three-phase Two-way Relaying 

Substituting the optimal rate assignment at each decoding 
region into the Lagrangian (l20l i. we have 



which gives 



n In / [l+73-PA(e)] ' P-).(7i,72,73) 



Pa. 

3 



X(i7id72d73 + / max { [l + 7iPA(e)] 

J-ies.2M-i 

[1 + -izPA{e) + 72Pfl(e) + 73/'A(e)72-Pi?(e)] 
X (71 , 72 , 73 )d7i '^72 rf73 

/ [l + 73^B(e)J ' P-y(7i,72,73) 

xdjidj2dj3 + / max { [1 +72 -PB(e)] ^ , 

[1 + IsPBie) + liPR{e) + 73/^3(6)71 /'fl (e)] } 
X (71, 72, 73)'^7i '^721^73^ + >^a{Pa - E^[PA(e)]) 

+\B(PE-^-y[PB{e)])+\R(P^-^^[PR{e)]), (60) 

where is the distribution function of 7. 

1) Region Mi (71 < 73,72 < 73)-' Setting the partial 
derivative at Pa{^) equal to zero, we can obtain 



^^-^[i+73P.(.r'^-A. 



TSKi I. "AVI 



^A73 
(701 



[l+73PA(e)] =Xa. (62) 



Hence, the optimal Pa{£) can be obtained as in ( l24l i. The 
optimal Pb{^) can be achieved using the same method. 

2) Region K2(7i > 73772 < 73)-' It is easy to see that the 
maximum Ra{^) is achieved when 

C(7iFA(e)) = C(73PA(e)) +C(72Pfl(e)), (63) 
which in turn gives 

(71 - 73)^^1 (e) 



Pni^) 



72 l+73/'A(e) 



(64) 



We substitute (|64] | into ( |60] |. take the partial derivative at 
Pa{^), and then obtain 



Ak(7i - 73) 



1 



-A, 



^2 [1+73^^(6)]' 

xP7(7ii 72, 73)^71^^72^73 = 0, (65) 

which yields ( l30l l. Therefore, the optimal PA(e) is given in 
( |29] l. The power allocation for Pb{^) is same as Ri. 

3) Region K3(7i < 73,72 > 73)-' We can use the sim- 
ilar method as M2 to get the optimal power allocation for 

PA{e),PB{e). 

4) Region R4(7i > 73, 72 > 73)-' According to (7) and (9), 
the maximum rates are achieved when 

C{jiPAie)) = C(73PA(e)) +C(72Ffl(e)), (66) 



C{j2PBie)) = C(73PB(e)) +C(7iPi?(e)), 
which offers 

(71 - 73)-PA(e) (72 - 73)-PB(e) 



Pni^) 



Therefore, 



72 [1 +73^14 (e)] 7i[l + 73/'s(e)] 



PAie) 



PBie) 



Pb 



r+(T-l)73PB(e)' 



rPAie) 



1 + (1-^)73^^(6)' 



(67) 



(68) 



(69) 



(70) 



(61) 



where t is given in (|33] |. If < r < 1, it always holds that 
if PA(e) > 0, then PB(e) > 0. Otherwise, if r > 1, then if 
Pb{^) > 0, it is sure that PA(e) > 0. In the following, we 
assume < r < 1, while the power allocation policies for 
r > 1 can be obtained in the same manner. 
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Substituting ( |68T l and ( fTOl i into ( |60] l and taking the the partial 
derivative of Pa{£), we can get 



, ojb Pb 



T72 



2 3 [l + (l-r)73PA(e: 

, _£b+3 

r72PA(e) 



1 + (1 - T)73P^(e) 



Ai?X7i - 73) 

72 



[1 + (1 - r)73PA(e)] 

;t fP7(7i:72,73)'^7ic'72'^73 = 0, 



[l + 73^A(e) 



(71) 



which can be simplified as equation i35[ . Thus, the optimal 
Pa is the solution of (|35T l as stated in ( [34l i. The optimal 
can be obtained similarly. 



Appendix B 
Derivation of Optimal Pa{£) and FB(e) for 
Two-phase Two-way Relaying 



According to the achievable rates in the decoding region 
ii,R2), we can rewrite (l42^ without loss of optimality as 



LOA 



-yew. 



liPA{e) 



(71. 72)^71 '^72 



l + 72Ps(e) 
-/ [l + 7i-PA(e)] ~P-^(7i, 72)^71^^72 ) 

UJB f f 1- — 

/ [l+72^s(e)J ' P-7 (71, 72)^71^^72 



72PB(e) 
l+7iP^(e) 

, 7i^A(e) 



^•7(717 72)^^71*^72 

2 

(71 7 72)^71*^72 



l + 72PB(e) 

Ha — — 

— [l + 7iP4(e)] ' ^-7(71, 72)^71^72 



fJ-B 

^J■B 



|_l+72PB(e)J '157(71772)^71*^72 



72Ps(e) 



, /3b 



Pt(71)72)c^7ic?7: 



l + 7iPA(e)_ 

+AA(Pl-E^[PA(e)]) +AB(P^-E^[PB(e)]), (72) 



by 



^A7l 



^A</)'i V7-T6Ri 2[l+72PB(e)] 

, 7i^A(e) ' ^ 



l + 72Ps(e) 

Pah 



P-y(71:72)c^7i<^72 



[l + 7iPA(e)] " p-y (71, 72)^710^75 
/3b7i 72 PB(e) 



■2 VAsRi 2[l+7iPA(e)]' 



72PB(e) 



l + 7iPA(e) 

Ma/3a7i 
2f?A[l + 72PB(e)] 



(71 > 72)0^71*^75 



1 + 72 Pi3(e) 

/iA/3yl7l 



-11a _ 



l + 7iPA(e) 



P-7 (71 : 72)0^71*^72 



/■ A^b/^b7i 72 PB(e) 
AeR^ 2(?B[l+7iPA(e)]' 

X P7 (71 > 72)0^71*^72 ) - Aa- 



PT(7i.72)oJ7iof72 

fin , 

72PB(e) 



l + 7iPA(e) 



(73) 



By differentiating on Pb{^) of (|72] |. the similar result can be 
obtained as 



t^A / f /^A7i72PA(e) 
i U-yeRi 2[l + 72PB(e)] 

7i^A(e) ' ^ 



l + 72Pi3(e; 

^B f f /3b72 



Pt(71:72)*^71*^72 



^2 VJtGR'i 

(71 1 72)0^710^72 



[l + 72Ps(e) 



/3s72 



7GR!, 2[l+7iPA(e)] 



1 



72PB(e) 
l + 7iPA(e) 

MA;3a 71 72 Pa (e) 



(71 7 72)0^71^72 
, 7i^'A(e 



heR[ 20A [1 + 72 PB(e)] 

XP7(7l:72)*^7l0^72 



l + 72PB(e) 

Mb/?-B72 

20b 



-1a _ 



x[l + 72PB(e)] - P7 (71 > 72)0^710^72 



Ms /3s 72 



2^5 [l + 7iPA(e)] 
(71 1 72)0^710^72 ) - A_B. 



72 PB(e) 



l+7iPA(e) 



(74) 



The partial derivative of (l72t with respect to PA(e) is given Let the derivatives equal to zero. If 7 G M.[, from ( f73T l. the 
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optimal condition should satisfy 



/ ujaPa _ ^MAA\ 71 



3 4 +2 

l+72^B(e) 



(75) 



which gives 



l + l2PB{e) 



l + l2PB{e) 



SXa- (76) 



From ( 1741 1. we can find the following optimality condition 



76R'i 



26 



B(p2 



26 f 



72 [l + 72-PB(e)] 



\29A(b[ 26 A ) [1 
, 7i^A(e) 



1 + 72^'B(e) 



71 72 -Pa (e) 
+ 72PB(e)]' 

P7(7l:72)'^7l'^72 0, (77) 



which offers 

72(a;B(/)2 ^ -MB)[l + 72-PB(e)] 
A*A)7i72-PA(e) 



[l + 72PB(e)]' 



7i^A(e) 



1 +72^*3 (e) 



= 0. 
(78) 



Replacing ( f76] l into ( fTsT l, and after some calculations, we can 
obtain the optimal power allocation of Pa{^) and Pb{^) in 
and 



Appendix C 

Derivation of Optimal Partition Criterion for 
Two-phase Two-way Relaying 



Here we only focus on the derivation for the threshold ^t^, 
while 7*'' can be obtained by the same way. Let we write the 
optimal 72^' as a function of 71, i.e., = /*(7i), where 
/*(7i) is the optimal function. We define /(71) = ./*(7i) + 
sg(7i), where s is any constant and 3(71) represents arbitrary 



variation. Thus, (l42l i can be rewritten as 

:^(/(7i)) 



--— in 



xd'yid'y2 



^B , 
— - — In 



00 /■/(71) 



Jo 

00 /'OO 



1 + 



liPA{e) 



^7(71,72) 



B \Ja Jo 



"'/(71) 
00 /•/(71) 



l+72^'s(e). 

1 -— \ 
l+7i-P4(e)J ^ P7(7ii72)c^7ic^72 1 



00 nOO 



"'/(71) 
00 /•/(71) 



MA 



^0 

00 poo 



l + 72-PB(e)J ^ P-y (71:72)^71 '^72 

fig 

P7(71j72)<^71C?72 

-9a 



72-PB(e) 
l+7iPA(e 



7i^A(e) 



1 + 72-PB(e) 



P-ihi,l2)d-iid'^2 



^^A 



[l+7iPA(e)] p-y (71, 72)^71^^72 

"'/(71) '^'4 

— l + 72PB(e) ' p-y (71, 72)^71^72 



JO 
00 />oo 



MB 



"'/(71) 



1 



72 -Ps(e) 



(71, 72)^71 '^72 

(79) 



l + 7iPA(e) 
+XA(Pl-^-y[PAie)]) + Ab(P^- lE^[PB(e)]) 

Intuitively, J'(/(7i)) gains its optimal value when s ~ 0, 
namely, the derivative of J'(/(7i)) is equal to zero when .s ~ 
0. Therefore, to obtain the optimal condition, it is necessary 
to satisfy fSO) 



s=0 



Then, it follows that 

UJA IJ-A 



-[l+7iP^(e)]- 



1 OA 

Ha ' 



0. 



7i^A(e) 



(80) 



l + ./*(7i)^B(e) 



.Pa. 



W(P2 OB 



[i + r{ii)PB{e)Y - 

xP7(7i,/*(7i))5(7i)rf7i = 0. 



r(7i)^B(e) 



l + 7iP4(e)J 



(81) 



Since the above equation needs to be hold for any 5(71), we 
can obtain 



1 



7lP-4(e) 
1+/* (71)^^23(6) 



\ - - - [l + 7iPA(e)]" 



.»B. 



6'a(wb02~^ - 11b) 



(82) 



6b{uja4>'i - fJ-A) 

Obviously, K > for all cases. Specifically, 7I'' = /*(7i) 
should be positive, otherwise it would not be well-defined in 
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